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We show that very small-amplitude oscillations of a highly symmetric, spheric or cylindrical, 
interface (thin membrane) between two fluids can result in inhomogeneous instability and breaking 
of the interface symmetry: the frequency of the breathing vibration selects the spatial symmetry. 
This mechanism may govern morphogenesis. 
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The nature of spontaneous symmetry breaking remains one of the most enigmatic questions of modern science. This 
problem emerges in connection with the equilibrium phase transitions, self-organization in noncquilibrium systems 
and many other areas in physics, chemistry and biology (see, e.g., Q), as well as with cell fission and morphogenesis, 
i.e., the development and spatial differentiation of complex body structures during their growth j2|. 

In 1952 Turing showed that the homogeneous state of some specific chemical reactions can lose stability with regard 
to a spontaneous increase of perturbations of certain form || . Since then the chemical basis is the prevalent idea of 
phenomenological theory of morphogenesis (see, e.g., Jl],||]). Turing's model is based on chemical or biological processes 
of local self-reproduction of some chemical agent (the activator) and far-ranging inhibition. As a consequence of such 
processes a very small increase of the activator concentration in a local region results in a global redistribution of the 
. substance concentrations and formation of more complex structure Q,^]. However, the Turing's chemical reactions 
are uncommon, unique and very complex processes. 

In this work we develop a new mechanism, without complexity, that breaks the symmetry by creating an instability 
O .■ in an oscillating interface, thin membrane, separating two different fluids. In other words, we show that if, for example, 
a spherical or cylindrical structure vibrates with a breathing symmetric mode for a given set of the frequencies the 
symmetry breaks with respect to bimodal, trimodal, pentagonal, etc. modes, i.e., the vibration frequency selects the 
, spatial symmetry of the interface. 
t-H ' We consider a thin symmetric membrane, spherical or cylindrical interface, with the radius Ro which separates 
two fluids with densities piand p2 {pi — P2 — p) respectively. Owing to Archimed's force the effective gravity 
acceleration operating on the internal fluid is g — g e (l — p\j p-i) « g e - We propose R is small enough, so the 
condition jk^/ p — jm 2 /RqP » g is valid. Here 7 is the surface tension and fc = m/R is the typical wave vector of 
the increasing deformation of the symmetric interface, m = 1, 2, 3, .... This is the condition when we can neglect the 
I gravity and consider only the effect due to the surface tension of the interface. 

Let us take, at first, for definiteness, a spherically symmetric interface S whose radius, R, oscillates with a frequency 
uj : R — Ro — dcosujt. From the incompressibility of the fluid it follows that its radial velocity is tvo = v m{t)R 2 \/ r ' 2 
where Vno(t) — dR/dt = dio smut. (This means that some source, for example, a small pulsating ball has to be inside 
the interface.) The vortex-free motion of an ideal liquid (we consider the effect of the viscosity below) is described by 
• the Euler and the continuity equations: 



dv dv , _. 1 

v 2 $ = v^ + vi$ = o (2) 

where $ is the velocity potential, v = and is the part of Laplacian depending only on coordinates of the surface 
S. For the undistorted spherical surface, from the symmetry of the problem, it follows that = 0, i.e., V^"^ = 0. 
Then, from Eq.|2|, we can write that = V r (w r o) = dv r o/dr + 2v r o/r = 0, in accord with tvo = «i?o(0^o/ r2 - I R 
the presence of a distortion, of the spherical surface S the interface radial velocity is v r = Vua{t) + ds/dt. Using 
this, we find from Eq.|l| that near the interface 

Z-W+*/"--%. (3) 
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where F(t) = du 2 cosujt is the acceleration of the interface and we neglect the term (viV)vi in Eq. |] by virtue of 
smallness of q [Q . Owing to smallness of q we can write the pressure near the surface as 

p = pF(t)(r-R-q) +1 \7lq+p (t) (5) 

Here we took into account that the pressure at the interface (when r = R + <r) is p — 7(171 + 02) +p (t) where a\ and 
(72 are the principal curvatures of the interface jjj: (01 + a 2) — Vj_<; since V^<r > Rq 1 - Substituting Eqj5|into Eq.[| 
we obtain 

^ = F(t)VxC - ^Vi? or |vi$ = F(t)Vi? - ^Vi? (6) 

We will seek solutions of the problem in the following form 

00 00 

M*)Sm and $ = J2 c m (t)t> m {r)S m - v m {t)Rl/r (7) 

m— m— 

where S TO is the complete orthogonal set of eigenfunctions depending only on the coordinates of the undisturbed 
surface S and satisfying the following equation 

(Vi + k 2 m )S m = (8) 

for r = Ro and the boundary conditions corresponding to the symmetry of the problem. In the spherical case S m = 
Ci,m-P)' m (cos#) exp(im<p) are the spherical functions of angles ip and 9 and — 1(1 + 1)Rq 2 where m = I, I — 1, — I 
and I = 0,1,2,... Substituting <£> from Eq.^ into Eq.|], using EqJ| and the condition V r (v r0 ) = cited above, we 
obtain the equation for 4 r m (r) : 

(V 2 . - k 2 J* m (r) = (9) 

with the boundary conditions V r ^ m — > when r — > and 4 , m (r) — A at r — Ro where A is some constant which 
does not reveal itself in the final results. Near the interface V r <& = v r = vrq(£) + dqjdt and so from Eq.f?] it follows 
that c m (t) — da m / dt(^I m )~^ Ro . Substituting $ from Eq.Q into EqJ^ and using Eq.y and c m (t), we find that 

00 

V i $ = - £ k 2 m >c m 1 S m da m /dt (10) 

m=0 

where x m = [V r ^' m /4' m (r)] r= ij; does not depend on the constant A . Then from Eq.[[o| and EqJ|, we obtain 

d 2 a m /dt 2 + [yk^xmp^ 1 - x m F(t)]a m = 0. (11) 
Using T = ivt/2 we can rewrite Eq.[n] as 

d 2 a m /dT 2 + (p m - 2q m cosu;t)a m = 0, (12) 

where 

q rn = 2x m d and p m = ti 2 ,^ 2 where fi^ = 4k 2 n ?c m 'yp~ 1 . (13) 

For the spherical interface k m Ro > 1 and x m ~ fc m = [Z(2 + l)] 1 / 2 i?Q 1 and so = 4[i(Z + l)] 3 / 2 i?Q 3 7P _1 and 
(Zm = 2d[Z(/ + l)] 1 /2 i? -i. 

These results can be extended easily to other cases. For example, when the interface have a form of a cylinder with 
vibrating radius, then S m = cos(fc;z) exp(imtp) and in Eq.[l3]>f m ~ k m and k 2 n — m 2 /Rq + ir 2 l 2 jh\ where ho is height 
of the cylinder. This vibrating cylindrical body can spontaneously distort in the axis z or with respect the azimuthal 
perturbations. 

We emphasize that Eq.|l2| coincides with Eq.(2.12) of Ref. || to describe the Faraday's instability |6) of the plane free 
surface of an ideal liquid under vertical periodic vibrations. These equations differ in the values of the parameters p m 
and q m . Moreover, in contrast to the Faraday's instability when the vibrations are reduced to trivial renormalization 
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of the gravity, in this work we consider spherical or cylindrical oscillating interfaces when the vertical direction, axial 
gravity, is not distinguished from other directions. Benjamin and Ursell S have constructed the stability diagram 
for Eq.[l2] with respect to the universal parameters p m and q m using the analogy between Eq.|l2| and the Mathieu's 
equation Q]. From this diagram it follows that the instability is realized only in regions near the points p m = n 2 
where n = 1, 2, 3, 4, .... In other words, the condition 

OJ = Lo n , m ~ n _1 f2 m = 2n~ 1 k m (x m ~f/p) 1/2 (14) 

determines the resonant vibration frequencies when the symmetric interface spontaneously deforms with respect to 
the standing wave with the azimuthal number m. However, the greater is n, the narrower is the width Eg (q m ) of the 
n-th region of the instability for given q m J^J^]. For the widest instability region, with n = 1, the value E g (q m ) ~ 2q m 
for q m < 1. It means that the instability takes place when (1 — q m ) < p m < (1 + q m ), he., the symmetry breaking is 
realized for the vibration frequency lying within the spectral range: 

fi m (l - x m d) < lu < Q m (l - x m d). (15) 



The threshold of the vibration amplitude d is limited by the fluid viscosity. For real fluid Eqsjll|and |12| include the 
additional terms j m da m (t)/dt and T m da m (T)/dT, respectively, where 7 m = 2vk 2 n C\ m and r m = Aisk^Cim/uj are 
proportional to the kinematic viscosity v and C m is some constant of the order of unity J§| . The threshold vibration 
amplitude, d — dt, for the instability region with n — 1, can be estimated from the condition E^\q m ) > 2T 

mi 1-6., 



q m = 2>c m d > T m . This condition follows practically from results of Refs. || and is obtained in [|9|. Using Eq.|15|, this 
condition can be written as 

d> d t = 2vC m k 2 n / x m uj ~ u(p/jx m ) ~ v(pR Q /^m). (16) 

For parameters of water dt — 4/im, i.e., the threshold vibration is a very small flutter of the interface. 

We propose that the results above may be used as a basis for a simple, without complexity, mechanism to trigger 
the fanciful morphogenesis appearing in nature. The frequency of homogeneous interface vibrations self-selects the 
space symmetry. If the interface oscillates with a characteristic frequency the germ symmetry will break when its 
radius Ro amounts to the quantity satisfied by Eq.[l5| or Eq.[l4]for n > 1. After the new symmetry appears the growth 
rate increases with surface curvature as is usual for many of Stefan- like problems [ flO| . 

The mathematical results reported here will be applied in a forthcoming paper to explain the morphogesis of the 
acetabularia, equinoderms and cell fision. 
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